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Abstract. Intermediate dispersion spectroscopic obser- 
vations of 37 Herbig Ae/Be stars reveal that the equivalent 
widths of their [OIJA6300 and Ha emission lines, are re- 
lated to their near-infrared colours in the same fashion as 
the T-Tauri stars. Such a correlation strongly supports the 
idea that the winds from Herbig Ae/Be stars arise in the 
same manner as those from T-Tauri stars, i.e. through ac- 
cretion driven mass-loss. We also find that the [OIJA6300 
line luminosity correlates better with excess infrared lu- 
minosity than with stellar luminosities, again supporting 
the idea that Herbig Ae/Be winds are accretion driven. 
If one includes the lower mass analogues of the Herbig 
Ae/Be stars with forbidden line emission, i.e. the classical 
T-Tauri stars, the correlation between mass-loss rate and 
infrared excess spans 5 orders of magnitude in luminos- 
ity and a range of masses from O.SM© to approximately 
IOMq. Our observations therefore extend the findings of 
Cohen et al. (1989) and Cabrit et al. (1990) for low mass 
young stars and, taken in conjunction with other evidence 
(Corcoran & Ray 1997a), strongly support the presence of 
circumstellar disks around intermediate mass stars with 
forbidden line emission. An implication of our findings is 
that the same outflow model must be applicable to these 
Herbig Ae/Be stars and the classical T Tauri stars. 
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1. Introduction 

Herbig Ae/Be stars (HAEBES) were first identified by 
Herbig (1960) as possible higher mass analogues of the 
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classical T-Tauri stars (CTTS). Subsequent work has 
shown that the two classes of young stars share many 
observational characteristics in common over and above 
those one might expect on the basis of observational selec- 
tion criteria alone. For example, like CTTS, the emission 
from a number of HAEBES is intrinsically and variably 
polarized (Jain & Bhatt 1995), many have IR excesses 
(Hillenbrand et al. 1992), water masers (Palla & Prusti 

1993) and evidence for accretion as seen in the redshifted 
absorption components of their NaD lines (Sorelli et al. 
1996). Moreover most HAEBES possess the signatures of 
strong mass loss, as indicated by P-Cygni line profiles 
(Finkenzeller & Mundt 1984; Perez et al. 1992; Imhoff 

1994) , broad forbidden emission lines (Bohm & Catala 
1994; Corcoran & Ray 1997a) and by their thermal radio 
continuum emission (Skinner et al. 1993). Estimates of 
their mass loss rates vary, and a number of different tech- 
niques have been used (for details see Nisini et al. 1995), 
but typical values are 10~* to 10~^ Mo/yr. 

The high level of activity of this group (for a general 
review of HAEBES activity the reader is referred to Catala 
1989) is perhaps initially surprising if one attempts to un- 
derstand it from a purely stellar perspective. If its source is 
the star alone then one would expect the mass loss to be 
driven by a sub-photospheric convection zone. However, 
according to conventional theory (e.g. Gilliland 1986) en- 
ergy transport in HAEBES is purely radiative. How then 
can one solve what Bohm & Catala (1995) call the "para- 
dox" of HAEBES activity? Two approaches have been 
taken (Bohm & Catala 1995): the first assumes that there 
is an inherent source of energy, such as internal rotation 
(Vigneron et al. 1990), which powers a stellar dynamo, 
while the second attributes the activity to viscous disk 
accretion in a manner similar to the CTTS (Hillenbrand 
et al. 1992). The viscous accretion disk model might, at 
first sight, seem more plausible given the observed IR ex- 
cesses of most HAEBES (e.g. Hillenbrand et al. 1992), 
however, alternative explanations for the IR excesses have 
been put forward, based on the idea that HAEBES are sur- 
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rounded by quasi-spherical clouds of gas and dust (Berrilli 

et al. 1992; Hartmann ct al. 1993; Mirosliniclicnko ct al. 
1997). At the same time, however, a number of HAEBES 
are associated with collimated optical and molecular out- 
flows (Mundt & Ray 1994) and given, what appears to 
be, the intimate connection between disks and outflows in 
CTTS (see, for example, Edwards et al. 1993) this would 
lead one to think that disks are present at least around a 
significant fraction of the HAEBES population. Thus the 
question still remains whether the standard disk model, 
developed to explain many of the observational features 
of CTTS (e.g. Bertout et al. 1988), is applicable to the 
higher mass HAEBES, or at least to those which, like the 
CTTS, have forbidden line emission (see also Ghandour 
et al. 1994; Corcoran & Ray 1997a)? 

If accretion disks, perhaps in association with stellar 
magnetic fields, provide the energy for driving HAEBES 
winds, one would expect a correlation to exist between the 
accretion rate, Maco and the wind mass- loss rate, Mwind- 
Such a correlation is found for the CTTS (Cabrit ct al. 
1990; Edwards et al. 1993), however, there are a number 
of important points that must be borne in mind before 
one can attempt a similar study for the HAEBES. Both 
Mace and M^ind are difficult to quantify and so appropri- 
ate substitutes have to be found. In the case of the CTTS, 
Cabrit ct al. (1990) and Edwards et al. (1993) used CTTS 
forbidden [01] line emission as a measure of wind strength 
and their infrared excess as a measure of accretion rate. 
Such an approach for CTTS can be justified since there is 
clear evidence that the forbidden [01] line emission orig- 
inates in a wind (e.g. Edwards et al. 1993) and because 
the infrared excess in CTTS correlates with the veiling 
continuum flux (see, for example, Hartigan et al. 1995), 
supporting the notion that the excess arises from viscous 
accretion. Can one use an identical strategy for HAEBES? 

In the literature there has been some question as to 
whether the forbidden line emission in HAEBES arises in 
the same manner as in CTTS (see, for example, Bohm 
& Catala 1994). A recent study by us (Corcoran & Ray 
1997a) however shows clear parallels between the forbid- 
den line emission as observed in HAEBES and CTTS (sec, 
for example, Hartigan et al. 1995). In particular double 
line profiles, consisting of high and low velocity blueshifted 
components, arc foimd in a number of HAEBES although, 
like CTTS (Hartigan et al. 1995), single low velocity 
blueshifted emission is by far the norm. The observation 
that the forbidden line emission of HAEBES is usually 
blueshifted, as with the CTTS, strengthens the conclu- 
sion that, it is formed in a stellar wind. Moreover the 
lack of any corresponding redshifted components also ar- 
gues that an obscuring disk is present (Corcoran & Ray 
1997a). We will therefore take the approach in this paper 
that the strength of the forbidden line emission can be used 
to gauge the strength of the wind from a HAEBES. There 
is some question whether the same approach can be taken 
in the case of the hydrogen emission lines. For example. 



it was thought in the past that the Ha line could be used 

as a gauge of wind strength; current modeling of Ha line 
profiles would suggest, however, that it is probably a bet- 
ter measure of accretion (Calvet 1997). In any event, we 
study both [01] forbidden and Ha line strengths here to 
see whether the same phenomenological relationships ex- 
ists between these quantities as in the CTTS. Moreover we 
also look at how both quantities depend on our more 'tra- 
ditional' measure of accretion, again to test the analogy 
with CTTS. 

Turning now to the measurement of Mace in HAEBES, 
this is somewhat problematic for various reasons. Accre- 
tion disk theory would suggest that we use either the ex- 
cess UV emission, which may arise from magnetic accre- 
tion columns close to the star, or that part of the infrared 
excess which derives from viscous dissipation in the cir- 
cumstcllar disk. The UV fluxes of these stars arc poorly 
known (although, indications of a far UV excess in a num- 
ber of cases is observed, Grady et al. 1993) so we are forced 
to consider the infrared excess. The infrared excess, how- 
ever, has two components which are difficult to disentan- 
gle: the first is due to reprocessed starlight in the disk 
and the second is due to heating by viscous accretion. For 
an optically thick, but geometrically thin, flat disk up to 
one quarter of the luminosity of the star, depending on 
the viewing angle, can be "reprocessed" (see, for exam- 
ple, Adams et al. 1987 and Strom et al. 1993). Moreover 
if the disk is "flared" (Kenyon & Hartmann 1987) even 
more light can appear as an "excess" although the con- 
tribution due to flaring is thought to be only important 
at mid-infrared or longer wavelengths. In the case of the 
CTTS samples (e.g. Hartigan ct al. 1995), one is dealing 
with a large number of stars of similar stellar photospheric 
luminosity but with widely different accretion rates, thus 
despite the "contaminating" effects of reprocessed light, 
statistically one can readily distinguish, for example, a 
correlation of the forbidden line strength with the accre- 
tion luminosity. The situation with HAEBES is not so 
clear-cut: HAEBES stellar luminosities can differ by fac- 
tors of 10^ (Hillenbrand et al. 1992) and one is dealing 
with statistical samples that are comparable in number 
to the CTTS samples. Thus a high degree of dependence 
on the stellar luminosity is introduced since, as one might 
expect, in absolute terms both the amount of reprocessed 
light and, to some extent, the accretion luminosity (see, 
e.g. Hillenbrand et al. 1992) tend to scale with the stellar 
luminosity. This tendency for both the accretion and re- 
processed contributions to increase with stellar luminosity 
in HAEBES, coupled with their large range in luminosi- 
ties, gives rise to, as we shall see, unavoidable correlations 
of the strength of the forbidden line emission not only with 
the infrared excess, as per CTTS, but with the stellar lu- 
minosity as well. A similar result was found by Nisini et al. 
(1995) when they studied a small sample of 14 HAEBES, 
although in their case they erroneously (see above) used 
the intensity of near-infrared HI recombination lines as 
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a measure of the wind strength. Obviously any correla- 
tion between the strength of the forbidden line emission 
and the stellar luminosity, can be used to support the 
idea, expressed by several authors (e.g. Bohm & Catala 
1994), that HAEBES winds are stellar rather than accre- 
tion driven. In the absence of a large sample of HAEBES 
within a narrow stellar luminosity range, statistically the 
question then becomes one of whether or not forbidden 
line emission is more tightly correlated with infrared ex- 
cess than with stellar luminosity. 

An alternative approach to circumvent any strong de- 
pendence on stellar luminosity is, in a sense, to "factor it 
out" as much as possible. In particular one can compare 
the equivalent width of forbidden line emission, and for 
that matter Balmer emission, with infrared colours since 
it is well known that the relative strength of the forbidden 
line emission, for example, increases with increasing red- 
ness in CTTS samples (Edwards et al. 1993). Moreover, 
Hartigan et al. (1995) found a clear correlation between 
the "veiling" index and K-L colour in CTTS, implying 
essentially a link between the ratio of the accretion lu- 
minosity and the stellar luminosity, Lacc/L*, and near- 
infrared colour. If the equivalent width of the forbidden 
line emission in HAEBES is then found to vary in the 
same way with colour as with the CTTS, then this would 
clearly point to accretion as the driving force for HAEBES 
winds. We initially adopt this promising approach here 
before considering any correlations of our wind measures 
with stellar or excess infrared luminosities. 

Finally, a relationship is also looked for between the 
strength of the HAEBES winds, as measured by the [01] 
line luminosity, and the amount of circumstellar dust 
present, as indicated by mass estimates from mm con- 
tinuum measurements. We find, in contradiction to the 
result of Nisini et al. (1995), who used a smaller sample, 
that there may be a weak correlation between the lumi- 
nosity of the [01] line and the amount of circumstellar 
matter. As will be described in the text, however, one has 
to be cautious in interpreting this result. 

After outlining the details of our observations and data 
reduction techniques in §2, we present the results of our 
study in §3 and discuss its consequences in §4. 

2. Observations and data reduction 

Our data derive primarily from three observing runs in La 
Palma on the Isaac Newton Telescope (INT) with the In- 
termediate Dispersion Spectrograph (IDS) from 28 August 
- 4 September 1991, 6-12 July 1993 and 14-20 December 
1994. Further observations were taken in La Silla with the 
ESO/MPI 2.2m Telescope using EFOSC-II from 23 De- 
cember 1991 - 1 January 1992. The instrumental details 
are listed in Table 1. 

The combined dataset from the three La Palma ob- 
serving runs and the La Silla run form a database of Ha 
and/or [OI]A6300 observations of 57 Herbig Ae/Be stars. 



Spec. 


Grating+CCD 


Disp. 


Res. 


Range 






(A pixcl^-"-) 


(kms"-') 


(A) 


IDS 


R632V+GE06 


0.70 


66 


400 




R632V+EEV5 


0.70 


66 


715 




R1200Y+GEC6 


0.36 


34 


205 




R1200Y-I-EEV5 


0.36 


34 


370 




R1200Y+TEK3 


0.39 


39 


500 




H1800V+GEC6 


0.22 


21 


120 


EFOSC 2 


Grism 9 










+Thompson CCD 


1.11 


105 


1145 



Table 1. IDS and EFOSC 2 grating-f-CCD parameters. The 
resolution (Res.) of the gratings and grism are calculated as- 
suming 2 pixels for full sampling. The pixel size of the Thomp- 
son CCD is 19/im, all other CCDs have a pixel size of 22/im. 
Data taken from ING User Manual VII and EFOSC 2 Oper- 
ating Manual. 



A subset was selected of this database of stars for which 
full photometric data was available in the literature (from 
the R to the N bands) and for which we had an exposure 
including the [OI]A6300 region. A total of 37 stars met 
these two criteria and are listed in Tables 2 and 3. The ex- 
cluded stars include both those for which we were unable 
to obtain (due to time constraints and/or weather prob- 
lems) an exposure of the [OI]A6300 region and those with 
photometric coverage insufficient to provide the necessary 
data for the calculation of the near infrared luminosity 
and excess. j-. 

Data reduction was carried out using standard IRAFEJ 
routines. Bias subtraction and flat-fielding corrections 
were determined from zero second exposures and tungsten 
lamp exposures respectively. The adjacent sky spectrum 
was subtracted from the object spectra, and the dispersion 
solutions were determined from CuAr arc exposures. 

Photometric data for HAEBES in the infrared and mil- 
limeter wavelengths were in most cases taken from the 
literature (Hillenbrand et al. 1992; Hamann & Persson 
1992) although it was found necessary to supplement these 
data with additional mm observations in the case of some 
stars. These millimeter continuum data were obtained at 
the Caltech Submillimeter Observatory (CSO) on Mauna 
Kea in April 1992. A ■^He-cooled germanium bolometer 
served as our detector with a beam size of 30" (FWHM). 
Uranus and 3C 273 were used as calibration standards. 
The CSO data was reduced as described in Beckwith & 
Sargent (1991). 

The [01] and Ha line luminosities were calculated 
from the absolute line equivalent widths and the contin- 
uum luminosity in the R band, centred on 7000A, using 
47rci^FAWiine, where d is the distance to the source. Fa 
is the flux density of the R band corrected for extinction 
and Wiinc is the absolute emission line equivalent width in 



^ The IRAF software is distributed by the National Opti- 
cal Astronomy Observatories under contract with the National 
Science Foundation 
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PV Ccp 


F2 


500 


7.0 


II 


13.80 


87.20 


1.15 


1.78 


0.52 


2.8 




Z CMa 


F5 


1150 


2.8 


II 


2.54 


24.80 


0.80 


1.66 


0.96 


1.69 




V645 Cyg 


AO 


3500 


3.1 


II 


6.00 


57.40 


1.98 


1.95 


1.25 


1.88 




BD+65 1637t 


B9 


1000 


1.9 


III 


<0.04 


12.71 


0.03 


0.51 


0.41 






LkHa 134 


B2* 


700 


2.1 




0.50 


31.54 


0.85 


0.70 


2.35 


2.57 





Table 2. The optical and IR data for the Herbig Ae/Be star sample. W([OI]) and W(Hq) are the absolute values of the equivalent 
widths of the [OI]A6300 and Ha lines respectively. The reddening corrected IR colours arc calculated from photometry published 
by Hillenbrand et al. (1992) and Hamann & Persson (1992). The classification by group (Grp.) is also from Hillenbrand et al. 
(1992). Those stars marked with (f) have only upper limits for W([OI]). The spectral type of those stars marked with an 
asterisk (*) is our own estimate, based on identification of absorption lines, where this differs from the value in the literature, 
see Corcoran & Ray (1996). All other spectral types, distance estimates and visual extinctions are from Hillenbrand et al. (1992) 
and Hamann & Persson (1992). Ha equivalent widths are with respect to the continuum level and not the expected photospheric 
absorption level. 



angstroms. This rncithod was used as a number of the ob- 
servations lacked suitable flux calibration standards that 
would allow a direct measurement of the line luminosi- 
ties. In the case of three stars (BD+40° 4124, BD+41° 
3731 and LkHa 134) no R band observations were avail- 
able and the V magnitude was used, extrapolated to the R 
band, e.g. from the SEDs in Hillenbrand et al. (1992). Es- 
timated errors for the luminosities are typically no better 



than ±20% depending on the uncertainties in the visual 
extinction and actual R band magnitude of the source at 
the time of observation. Uncertainties in the distance and 
visual extinction estimates can introduce sizable errors, 
see Discussion. Estimates of the visual extinction, partic- 
ularly, vary greatly in the literature (e.g. LkHa 198 with 
a quoted Av of 2.5 in Hillenbrand et al. 1992 and an Av 
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Fig. 1. Absolute equivalent widths (in angstroms) of [OI]A6300 and Ha plotted against near-infrared colours, H-K and K-L. In 
this and subsequent figures, Hillenbrand et al. (1992) Group 1 stars arc indicated by diamonds. Group H stars by crosses and 
Group HI stars by squares. Note that the Group HI stars cluster around the expected IR colours of main sequence stars, while 
the Group I and II stars all possess H-K and K-L colours equal or greater than those expected from an opaque geometrically 
thin circumstellar disk. The arrows indicate that the value of L[oi] is an upper limit only, and [OI]A6300 was not detected from 
that star. A K-L colour is not known for the Group I star LkHa 218 and so it is missing from the plots to the right. Note also 
that one Group HI star, BD-l-41 3731 was seen in absorption in Ha and is therefore not shown in the bottom plots. It is almost 
certainly not a pre-main sequence star (see text) 



of 5 in Berrilli et al. 1992) and. reflect the uncertainty in 
the dust composition and distribution about these stars. 

Table 2 lists the 37 HAEBES selected from our dataset 
according to the criteria described above. Of these 24 
show [OIJA6300 emission. The remaining 13 stars have no 
[OI]A6300 emission and are marked with a (f) and the de- 
tection limit quoted is an upper limit for W([OI]). These 
will be discussed later. The broad band photometry val- 
ues in Table 2 are taken from Hillenbrand et al. (1992) 
with the exception of the stars PV Ccp, Z CMa, V645 Cyg 
and LkHa 134 which are from Hamann & Persson (1992). 
All are corrected for interstellar extinction, using the vi- 
sual extinctions quoted in Hillenbrand et al. (1992) and 



Hamann & Persson (1992), assuming the standard inter- 
stellar reddening curve with Ay/Eb-v ~ 3-1, and that 
the ratio of the absorption values. Ax /Ay scale, with 
wavelength, as tabulated in Scheffler & Elsasser (1982). 
Table 3 presents the calculated luminosities of [OI]A6300 
and Ha, along with the infrared and stellar (photospheric) 
luminosities for the 37 HAEBES. The V and R magnitudes 
of the stars as well as the bolometric correction applied in 
each case is also listed. 



3. Results 
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BD+65 1637t 
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12.00 
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Table 3. Table of observed broadband V and R magnitudes plus bolometric correction (B.C.) applied in calculating the 
photospheric luminosity, L*, of individual stars (see text). The table also gives the derived line luminosities, Lire the infrared 
excess (over the range 0.7/jm to 10.2/Lim, see text for precise definitions), L,,nir the expected infrared contribution from the 
stellar photosphere (over the same range) as well as the photospheric luminosity, L,, for the Herbig Ae/Be star sample. All 
luminosities are corrected for reddening. Those stars marked with (f ) have only upper limits for L[oi] • ABS indicates the line is in 
absorption. The bolometric corrections are taken from Zombeck (1990) and the V and R magnitudes are taken from Hillenbrand 
et al. (1992) and Hamann & Persson (1992). 



3.1. Relative wind strength and near-IR colours 

A number of infrared colours were chosen to sample the 
dust temperature gradient and, indirectly, the accretion 
activity; H-K, K-L, L-M, M-N, and IRAS 12 and 25/im 
fluxes converted to magnitudes (see Table 2). The near- 
infrared photometry (H; 1.65/xm, K; 2.2/Km, L; 3.6/xm, M; 
4.8/um, N; 10.2/im) and the IRAS observations at 12/im 
and 25/im effectively sample different temperature regimes 



of the circumstellar dust, ranging from about 1700K (i.e. 
close to the dust sublimation point) at 1.65/im to about 
lOOK at 25/im (these are approximate values from Wein's 
displacement law) . The IR colours are accurate to approx- 
imately ± 0.02 mag. 

Fig. 1 plots the relationship between the absolute 
equivalent widths of [OIJA6300 and Ha and the near-IR 
colours H-K and K-L, given in Table 2. The HAEBES are 
clearly separated into three distinct regions of the plot de- 
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Fig. 2. a) Near infrared colour K-L plotted against the normalized infrared excess luminosity, Liro/L,. The infrared colour 
increases with increasing near-infrared excess luminosity, as expected from the connection found by Hartigan et al. (1995) 
between K-L colour and the veiling luminosity in CTTS. The Group I stars to the left of the main distribution are discussed in 
the text, b) Absolute equivalent widths of [OIJA6300 and Ha. The forbidden line scales well with the Balmer line. Those Group 
I stars without detectable [OIJA6300 (identified in Table 2) also show the weakest Ha equivalent widths of all Group I stars, 
suggesting reduced activity in the wind is measured by both the Ha and [OIJA6300 lines. The Group H stars appear to have 
the highest levels of [OI]A6300 for a given Ha which is in keeping with their less evolved state, more active circumstellar disks 
and consequently more powerful outfiows. 



pending on their Hillenbrand et al. (1992) group (see Table 
2). Here wc recall that Hillenbrand et al. (1992) classified 
HAEBES by the slope of their near-infrared spectral en- 
ergy distribution (SED). Group I stars show SEDs with 
a declining slope in the near-infrared that can be mod- 
eled with a geometrically flat optically thick disk with a T 
oc r~^/^ temperature law. The Group H stars have SEDs 
with flat or rising slopes in the near to mid-infrared, and 
the Group HI Herbig stars have near-infrared SEDs with 
little or no deviation from a normal stellar spectrum. The 
Group III stars in our sample all show values of H-K and 
K-L close to the values expected for main-sequence B-type 
stars (~ 0.0) and no detectable [OIJA6300 emission. Their 
distinct separation from the Group I & II stars in Fig. 
1 is analogous to the division between the weak-line T 
Tauri stars (WTTS) and the CTTS in the corresponding 
plots for lower mass stars (Edwards et al. 1993). All of 
the Group III stars in our sample show absolute equiva- 
lent widths of Ha < 15 A. The possibility of the Group HI 
stars being an intermediate mass analogue to the WTTS 
is discussed later. The cut-off colour indices for the Group 
I stars are H-K = 0.4 and K-L = 0.8 in line with the 
Hillenbrand et al. (1992) definition of this group (see also 
§4.1). The transition from Group I to Group II stars is 



then one of increasing redness as the relative contribution 
from the near-infrared to the SED goes up. 

A clear result in Figs, la and lb is that the equiv- 
alent width of [01] emission, W([OI]), scales with H-K 
and K-L and, in particular, the 3 groups have very dif- 
ferent [01] emission line properties. While the Group HI 
stars, i.e. those without optically thick disks (H-K < 0.4 
and K-L < 0.8), have no detectable [01] emission, the 
equivalent width of the [01] emission for Group I stars 
rises by over 2 orders of magnitude with increasing K- 
L. The Group II YSOs then constitute those stars with 
the largest [01] equivalents widths. The results for the 
W(Ha) vs. H-K & vs. K-L (Figs. Ic, d) are found to be 
comparable to those for W([OI]). A similar effect was dis- 
covered by Ghandour et al. (1994) for a smaller sample 
of HAEBES when looking for a correlation between the 
equivalent width of [OI]A6300 and I-N colour. As Ghan- 
dour et al. (1994) were comparing the ratio of the flux 
at 0.9/um (I) and 10.6/zm (N) they did not find as clear 
a separation (for reasons that will be discussed shortly) 
between the Group III stars and the Group I and II stars 
as we do here. 

A very interesting finding, if we compare Figs, lb and 
Id with Fig. 2 of Edwards et al. (1993), is that the de- 
pendence of both [01] and Ha equivalent width on K-L is 
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Fig. 3. a) Emission line luminosity of [OIJA6300 plotted against the infrared excess, Lire- b) Emission line luminosity of [OI] 
A6300 plotted against photospheric luminosity, L* . c) Emission line luminosity of Ha plotted against the infrared excess (Lire)- 
d) Emission line luminosity of Ha plotted against the photospheric luminosity, L«. All luminosities are in units of L©. See Fig. 
1 for key. 



virtually identical in both the HAEBES and CTTS. As in 
the CTTS, the absolute value of the [01] equivalent width 
is found to vary in the HAEBES by approximately two 
orders of magnitude, from 0.1 to lOA, as K-L is increased 
from 0.5 to 1.5 while the corresponding relative change 
in the Ha equivalent width, from 10 to about lOOA, is 
smaller, only one order of magnitude. Such close parallels 
in phenomenological relationships argue strongly that the 
same basic accretion/outflow process is operating in both 
groups of stars. 

The good degree of correlation between [O1JA6300 
equivalent width and near-IR colours for the HAEBES 
does not extend to colours longward of ^4.0/im. More- 
over in the corresponding plots at longer wavelengths, the 
3 Hillenbrand et al. (1992) groups do not occupy distinct 
regions of IR colours as in Fig. 1. For example, compar- 



isons of the W([OI]) vs. L-M show little correlation of the 
two variables. The circumstellar conditions giving rise to 
the L-M, M-N and IRAS 12/xm - 25/xm colours must de- 
part significantly from the conditions in the inner region 
(as indicated by the H-K and K-L colours). Note that the 
material emitting at ~ 4/Ltm is typically situated at about 
an AU from an AO star (assuming a flat optically thick 
reprocessing disk with an r~^/^ temperature law). [Table 
4, which will be discussed later, shows the results of two 
survival analysis tests on the W([OI]) vs. IR colours data. 
Both tests agree with the visual correlations seen in Fig. 1 
and confirm the lack of correlation between W([OI]) and 
L-M, M-N & IRAS 12-25/im colours.] 

Fig. 2a plots K-L colour against the normalized in- 
frared GXCGSS, Lire /L* . Lire L* are listed in Table 3 
and are defined more precisely below, but can be taken as 
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Fig. 4. a) The infrared excess luminosity, Lire vs. the photospheric luminosity, L«. b) Emission line luminosity of [OI]A6300 vs. 
Ha. All luminosities are in units of Lq. See Fig. 1 for key. 



the excess IR luminosity roughly over the range 1-10/im 
and the stellar bolometric luminosity, as determined using 
the reddening corrected visual magnitude, respectively. 
From Fig. 2a it can be seen that K-L colour is correlated 
with Liro/L* and that the correlation extends to values 
of Lirc/L* well beyond the expected regime for pure re- 
processing disks. Moreover an identical correlation (e.g., 
Edwards et al. 1993) with approximately the same slope 
and a similar range in (K-L, Ljrc/L*) space is found for the 
T Tauri stars (TTS). This clearly suggests that, as with 
the TTS, K-L is a measure of the normalized accretion 
luminosity (Lacc/L*)- 

In Fig. 2b we have plotted our relative measure of wind 
strength, W([OI]), against W(Ha). Not surprisingly the 
two equivalents widths appear correlated and the corre- 
lation is exactly the same as found with the CTTS (e.g., 
Cabrit et al. 1990) with the same typical ratio between the 
[01] and Ha line strengths. In particular Fig. 2b also shows 
that those Group I stars with no detected [OIJA6300, have 
low absolute Ha equivalent widths. 

3.2. Luminosity measurements 

For the purposes of the analysis here, a number of quan- 
tities are defined. The observed near-IR luminosity, Lnir, 
is the luminosity of a given star over the spectral range 
0.7/Ltm - 10.2/im. The corresponding photospheric lumi- 
nosity, nil-, is calculated as the integrated emission from 
a blackbody of temperature Teff find radius as given by 
Hillenbrand et al. (1992). Consequently, the infrared ex- 
cess, Lire, is = Lnir " L*,nir- This is the value we use as a 
diagnostic of the accretion luminosity. The stellar luminos- 



ity, L* quoted in Table 3 is calculated from the redden- 
ing corrected visual magnitude of each star, with visual 
extinctions from Hillenbrand et al. (1992) and Hamann 
& Persson (1992), using U/Lq = io-0-4(Mboi-4.72) 
main-sequence bolometric corrections (Zombeck 1990). 

All the Lnir luminosities were calculated, assuming 
spherically isotropic emission, by a trapezoidal integra- 
tion of the published photometry (Hillenbrand et al. 1992; 
Hamann & Persson 1992) over the spectral ranga covered 
by the bands R to N (roughly 0.7/im - 10.2/im) El. 

The results of the various luminosity correlation inves- 
tigations are plotted in Figs. 3, 4 and 5. Figs. 3a, b com- 
pare the [OIJA6300 emission line luminosity with Lire and 
L* respectively to determine the relationship between our 
diagnostic of wind mass-loss rate with the accretion and 
stellar luminosities respectively. Figs. 3c, d are the corre- 
sponding plots for the Ha line luminosity. Fig. 4 shows 
the Lire vs. L* relationship and the comparison of the two 
emission line luminosities against each other. Fig. 5 is a 
combined plot showing L^qi] line luminosity against Lire 
for our sample of HAEBES and the TTS sample of Cabrit 
et al. (1990). 



This wavelength range was chosen because any signature 
of accretion is more likely to be dominant in the near-infrared 
whereas one expects a proportionally increased contribution 
from reprocessed light at longer wavelengths. This should be 
particularly notable in the case of a flared disk (Kenyon & 
Hartmann 1987) where the contribution from reprocessed light 
at 60/im with a surface scaling as i?^ can be increased by a 
factor of 10, for z = 9/8 or a factor of 33, for z — 5/4, over 
that from a flat disk (z=0). 
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Fig. 5. Emission hno luminosity of [OI]A6300 plotted against the infrared excess (Lire) for our sample of Herbig Ae/Be stars 
and the T-Tauri stars from Cabrit et al. (1990). All luminosities are in units of L©. TTS are represented by x's, all other stars 
are represented as described in Fig. 1. Those TTS with arrows are cont/e sources as described in Cabrit et al. (1990) and may 
have incorrectly calculated visual extinctions. 



The data points in Fig. 3 arc separated by group (ac- 
cording to Hillenbrand et al. 1992), with downward arrows 
indicating those stars with upper limits to the observed 
[OIJA6300 emission. One can see that there is a clear re- 
lationship between L^qi] our wind tracer, Lua, and Lire 
(Figs. 3a and c). As discussed in the Introduction, we 
expect, and find, that L[oi] and Lho also correlate with 
L* although the correlations (Figs. 3b and d) are not as 
good as with Lire (see also the results of the statistical tests 
which are described below). Note, once again, that all of 
the Group III stars show no detected [OIJA6300 omission 
and have only upper limits plotted in Fig. 3. Fig. 4a shows 
the expected correlation between Lire and L* although we 
note that the Group II stars show the greatest scatter in 
the plot reflecting their wide range of accretion luminosi- 
ties for a given L*. In Fig. 4b we sec that L^qi] is typically 
1% of Lho underlying the point, already made by Fig. 2b, 
that, while we might expect some forbidden line emission 
even from those stars with weak Ha emission, it will be 
difiicult to detect. Fig. 5, shows a plot of L[oi] against Lire 
again but this time we have included the corresponding 
data for 36 TTS from Cabrit et al. (1990). Note that the 
Lire calculated here is equivalent to the near to mid-IR 
excess luminosity, L^ir as defined in Cabrit et al. (1990). 
Comparing the data for the two classes of stars we find 
a smooth continuation of the correlation between the two 
variables over a luminosity range of 5 orders of magni- 
tude. The mass range sampled covers ^ O.SM© to approx- 
imately 10 Mq. In Cabrit et al. (1990) and Fig. 5 there 



are four TTS marked with arrows indicating that they are 
continuum sources (cont/e) in Cabrit et al. (1990) which 
may have incorrectly calculated visual extinctions. Cabrit 
ct al. (1990) note that decreasing their Ay by 1 moves 
them in the direction indicated by the arrows, reducing 
the degree to which they deviate from the trend shown by 
the HAEBES and the TTS. 

The visual impression of correlations is supported by 
two statistical survival analysis tests. The Generalized 
Kendall's Tau and Cox's Proportional Hazard tests de- 
termine the probability that a random set of uncorrelated 
(x,y) points show the same observed degree of correla- 
tion (see, for example, Isobe et al. 1986). The tests deal 
with censored bivariate data and, therefore, the stars with 
no detected [OIJA6300 emission arc included in the set as 
having line emission equal to the upper limits determined 
from the minimum detection level in each case. The re- 
sults of the survival analysis tests are shown in Table 4, 
and indicate that the likelihood of the observed correla- 
tion occurring by chance for L[oi] and Lhq vs. Lire is very 
small. 

It is interesting to note that the tests show the Lua 
vs. Lire correlation and the Lua vs. L* correlation are 
statistically equally probable, whereas the L[oi] vs. Lire 
correlation is much tighter than the correlation between 
L[oi] and L*. We had already noted this to be the case on 
the basis of visual inspection. The fact that the relation- 
ship between the line luminosity ^[oi] and the stellar lumi- 
nosity is weaker than that between L[oi] and Lire clearly 
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0.0025 
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Table 4. The results of the statistical survival tests. All probabilities are those of no correlation existing between the compared 
variables. For comparison the probabilities of no correlation existing for the sample of 36 TTS in the sample of Cabrit et al. 
(1990) is < 0.0001 for the comparison of L[oi] and Lhq against both Lire and L*. The luminosity range sampled by Cabrit et 
al. (1990) was of course smaller. 



suggests that Mwind depends on Mace- In particular, the 
root of the larger scatter in the L[oi] vs. correlation 
lies in the huge increase in W[OI] by over two orders of 
magnitude with increasing K-L. This creates a separation 
between the 3 groups at a given L*, in the sense: L[oi] 
(Group II) > L[oi] (Group I) > L[oi] (Group III). Con- 
versely, the much smaller difference in tightness between 
Lhq vs. Lire and Ljjq vs. L* is a direct consequence of 
the much smaller span in W(Hq!) between the 3 groups, 
compared with W[OI]. 

3.3. Millimeter measurements 

A total of 16 Herbig Ae/Be stars with [01] measurements 
were observed with the CSO in 1992. Two of the stars, 
LkHa 257 and LkHa 259, lacked complete optical and 
near-infrared photometric data and are therefore not in 
Tables 2 and 3. Henning et al. (1994) made millimeter 
observations of an additional 2 stars for which we have 
[01] data. Table 5 gives the combined sample of 18 stars, 
with the relevant continuum flux at either 1.1 mm or 1.3 
mm depending on the observations. By assuming that any 
mm emission is optically thin, it is possible to estimate an 
upper limit for the mass of the circumstellar material, Mcs, 
(see, for example. Mannings 1994). 

Mcs « (1) 



Here F\ is the flux density at wavelength A, d the dis- 
tance to the source, Bx (T) the Planck function and kx is 
the mass absorption coefficient, kx is taken to be a power 

law characterized by kx — kq {^)^ ■ The value of kq is as- 
sumed to ^ 0.1 cm^g^^ at Xq = 0.25 mm, as used by Man- 
nings (1994), Beckwith et al. (1990) and Beckwith & Sar- 
gent (1991). Following Mannings (1994) we take the value 
of /? ~ 1.0 The masses quoted in Table 5 assume a gas to 
dust ratio of 100:1 (this factor is contained in the normal- 
ization of Kx) and are thus total gas+dust masses. A dust 
temperature of 37 K was assumed in all cases (Mannings 
1994). These mass estimates are subject to some varia- 
tion due to the differing beam sizes used in the different 
observations. Henning et al. (1994) used a 23" beam, as 
opposed to the 30" CSO beam. The difference in beam 
size will result in varying amounts of background contam- 
ination in the flux measurements and consequently over- 
estimates of the target's circumstellar mass. This effect is 
most pronounced for the most distant stars where the re- 
gion sampled is much larger than that for the nearer stars. 
For example, V645 Cyg is calculated as having a circum- 

^ Note however that Pollack et al. (1994) suggest that /3 ^ 1.5 
over the millimeter range. Combined with their values for kq, 
the resulting k\ would be a factor of approximately 2 smaller, 
resulting in somewhat larger estimates of the circumstellar 
masses. The calculated masses may need to be revised upward 
should the work of Pollack et al. (1994) prove more accurate. 
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Table 5. Millimeter continuum measurements and circumstcUar gas+dust mass estimates. The data is from CSO observations 
at 1.1mm, with the exception of the two stars marked with an asterisk (*) which are 1.3mm observations from Henning et al. 
(1994). The mass is calculated from Equation 4 and the error is of order 20%. There are two mass estimates for V645 Cyg (see 
text for details). 



stellar mass of nearly 9OM0, which is clearly an overesti- 
mate. To compensate for this we assumed a 1/r^ density 
law for the dust and corrected the 1.1mm emission from 
V645 Cyg to the average distance of the rest of the sample. 
Note that even with this compensating factor the cir- 
cumstellar mass of V645 Cyg is still anomalously high. At 
the distance of V645 Cyg (3.5 kpc) the 30" beam encom- 
passes nearly a 0.5 pc diameter. There may be multiple 
sources in this region (see Corcoran & Ray 1997b), each of 
which might contribute to the 1.3 mm flux. Table 5 shows 
the two values for the V645 Cyg circumstellar mass. The 
masses listed in Table 5 are best treated as very rough es- 
timates only although we note that for the many stars we 
have in common, our mass estimates agree with those of 
Mannings (1994). There are several reasons why we should 
treat the quoted values with suspicion: firstly as we have 
already stated opacity values at millimeter wavelengths 
arc at least uncertain by a factor of 10 and besides they 
may vary from place to place (Men'shchikov and Henning 
1997). Moreover assuming the dust distribution around a 
young stellar object (YSO) is isothermal is clearly artifi- 
cial, one should use much more sophisticated models (Di 
Francesco et al. 1997). 

Figs. 6a and b show plots of the observed [01] luminos- 
ity and Lire against derived circumstellar mass. Although 
the sample is relatively small, there appears to be a weak 
correlation between Lpi] and Mcs and Lire and Mcs- We 
will consider the possible origin of this correlation in our 
Discussion. 



4. Discussion 

Dust colours 

The throe groups of HAEBES, as defined by the shape of 
the spectral energy distributions (SEDs) (Hillenbrand et 
al. 1992) fall into distinct regions in the W([OI]), W(Ha) 
vs. H-K, K-L plots. The correlation observed between the 
relative strength of the forbidden emission and the noar-IR 
colours for Groups I and II indicate a possible relationship 
between the relative strength of the wind and the tem- 
perature gradient in the circumstellar material, with the 
equivalent width of [OIJA6300 in Group I varying by over 
2.5 orders of magnitude with increasing K-L. This cannot 
be easily understood in terms of purely stellar processes, 
and points to some influence of circumstellar environment 
on the formation of the line. A plausible explanation for 
this behaviour, as we shall now demonstrate, is the pres- 
ence of an optically thick disk (or a combination of disk 
plus halo) around Group I and H stars. 

The expected near-IR colours from an optically thick, 
passive reprocessing disk or flat accretion disk (Bertout et 
al. 1988) can be calculated simply from the assumption of 
a AFa oc A~^/^ relationship between the flux density and 
wavelength: 

- (r^)"'' (2) 

f\{mK) Xk 

where fxim^) 
is the spectral irradiance in erg cm~^s~^A~^ for an ob- 
ject of magnitude rux in the x fllter. Here Xh = 1.65/im 
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Fig. 6. Emission hnc luminosity of [OI]A6300 (a) and Lire (b) plotted vs. the calculated circumstellar mass of gas and dust, in 
units of Lq and M0 respectively. Symbols are as in Fig. 1. 



and \k = 2.2/im are the central wavelengths of the H 
and K filters respectively. Using the zero magnitude ir- 
radiance figures provided in Zombeck (1990), one finds a 
value of thh - ts^k = 0.4, or in commonly used notation, 
H-K = 0.4. Similar calculations for K-L yield a value of 
0.82. For most of the Group I stars, and all of the Group 
II stars, the values of H-K, K-L exceed those possible from 
a flat reprocessing disk and three obvious explanations are 
possible. Firstly the disk could be fiared in some manner 
increasing the amount of intercepted flux (see, for exam- 
ple, Kenyon & Hartmann 1987) and flattening or raising 
the 1.65-3.6 jim portion of the SED. Flaring of the disk 
seems to us to be an unlikely explanation because of ge- 
ometry: even in a very flared disk, the inner portion (r < 
10 AU; A « 1-10 /Ltm) should remain relatively flat in com- 
parison to the outer disk which can intercept significantly 
more stellar photons. Another possibility is the presence 
of an unresolved IR companion that distorts the true near- 
infrared SED of the star. While this might be important 
in a small number of cases, recent searches in the infrared 
for HAEBES binaries (e.g., Leinert et al. 1997; Pirzkal et 
al. 1997) show that any companions are usually faint in 
comparison to the HAEBES and so cannot explain their 
infrared excesses. The final possibility is that the increased 
reddening is due to a dusty halo which scatters starlight 
back onto the disk and heats it at large radii (Natta 1993; 
Mannings 1994; Di Francesco et al. 1997). If one accepts 
that YSOs with halos plus disks tend to accrete at higher 
rates than comparable stars with disks alone, then the halo 
plus disk model explains not only the findings of Hartigan 
et al. (1995), who showed a correlation between the in- 



frared colours and veiling in TTS but, by extension, our 
result that the relative strength of the forbidden line emis- 
sion is correlated with infrared colour in HAEBES in the 
same way. 

Looking at Figs, la and c, there are, however, a small 
number of stars (approximately half a dozen) that do not 
fit the general trend, i.e. they have very red colours but no 
detected forbidden line emission. Almost exclusively, they 
are Group I stars. In some cases, we may not be deal- 
ing with a bona fide HAEBES. For example, HD 163296 
(MWC 275) is not in a star forming region, it does not ap- 
pear to be associated with any reflection nebulosity, and 
its spectral characteristics are inconsistent with it being 
a HAEBES (The et al. 1985). Di Francesco et al. (1997) 
question whether MWC 137, originally miss-classified as 
a planetary nebula, is, in fact, a pre-main-scquence star. 
In some cases variability of the forbidden line emission 
may be the reason why it was not detected by us. For 
example, although we did not observe [OIJA6300 emission 
from LkHa 218, Bohm and Catala (1994) did at 3 times 
our detection level. Variability in forbidden line emission 
strength, by factors of 200-300% has also been found in a 
number of other HAEBES (Bohm and Catala 1994; Corco- 
ran and Ray 1997a). Irrespective of their infrared colours, 
Fig. 2b shows that the Group I stars with no observed 
forbidden line emission, are also weak Ha emitters. This 
may imply a link between the current accretion and out- 
flow rates: Hartmann et al. (1994) have shown that for 
TTS the bulk of the Ha line emission comes from infalling, 
magnetospheric material and is not formed in the wind. 
Thus the Ha equivalent width, as we mentioned in the 
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Introduction, is probably better regarded as a signature 
of accretion than of wind strength. This begs the ques- 
tion, could some of the Group I stars without detected 
[01] emission be undergoing episodes of reduced accretion 
activity? If this is the case then their near-infrared colours 
could be due to the presence of a optically thick reprocess- 
ing disk which would give H-Krs 0.5. 

Finally we remark in connection with the other group 
for which we did not detect forbidden line emission, i.e. 
the Hillenbrand Group 3 stars, that a significant fraction 
may also not be genuine HAEBES. For example HD 37490 
(Omega Orionis) and HD 76534 are almost certainly clas- 
sical Be stars (Adelman 1992; Oudmaijer and Drew 1997). 
Similarly BD -1-41 3731 was also rejected as a member of 
the HAEBES class (The et al. 1994) consistent with our 
observation of Ha in absorption in this star. 

In summary, the relationship between the [01] and Ha 
equivalent widths with near-IR colours appears identical 
to that found for the TTS. It suggests, as with the TTS, 
that while some HAEBES are surrounded by disks, others 
have a spectral energy distribution that derives from a 
combination of a disk plus a halo. In a number of cases, 
the Group III stars, either no disks are present or at best 
they are optically thin, and it is precisely these stars that 
show very little sign of accretion and outflow activity. In 
contrast, their much redder counterparts, the Group II 
stars show strong Ha and [01] emission. 

4-2. Luminosity relationships 

Figs. 3a and 3b show good correlations between Ljqi], Lhc« 
and the infrared excess. Lire- As with the TTS, this sug- 
gests a link between Mwind and Mace in the HAEBES. 
When combined with the results from Cabrit et al. (1990), 
in Fig. 5, the correlation between wind and accretion, in 
the form of the L[oij vs. Lire plot, spans 5 orders of mag- 
nitude in Lire- It is difficult to imagine a purely stellar 
physical process that could produce such a correlation 
especially when one considers that the sub-photospheric, 
photospheric and chromospheric conditions that occur in 
TTS are thought to be very different in HAEBES (for 
example, Catala 1989; Palla & Stabler 1993). It is partic- 
ularly interesting that the L[oi] vs. Lire plot (as well as the 
Lfja VS. Lire plot which is not reproduced here) smoothly 
connects the two stellar groups. This clearly suggests that 
disk plus halo models, which have met with much success 
in explaining many of the characteristics of CTTS, must 
be appropriate in a "scaled-up" version to their somewhat 
more massive counterparts. 

Other researchers' investigations of possible relation- 
ships between M^ind and Mace have produced similar re- 
sults. Hamann & Persson (1992) used the strongest line of 
the calcium IR triplet, CaIIA8542, as a measure of Mwind 
in both TTS and HAEBES and compared it with the 
infrared excess calculated over a wavelength range 1 - 
25/im, using ground based photometry and IRAS data. 



The CaIIA8542 luminosities scale with the infrared ex- 
cess in the same manner as the emission line data pre- 
sented here. Hillenbrand et al. (1992) showed a correla- 
tion between Lna and Maee for both TTS and HAEBES 
and found a result similar to Fig. 5. Hillenbrand et al. 
(1992) used the monochromatic IR excess at 3.5/zm, tak- 
ing into account the passive reprocessing of starlight by 
the disk material, as a measure of Mace- As pointed out 
earlier on, however, current models would imply that Lho 
is not a measure of wind strength and that any calcula- 
tions of mass loss rates based on the flux of the hydrogen 
lines, both optical and infrared, are "meaningless" (Cal- 
vet 1997). This criticism would also apply to the work 
of Nisini et al. (1995) who claimed to have directly mea- 
sure Mwind in HAEBE stars using infrared HI recombina- 
tion lines (Pa/3, Br7, Bra and Pf/3). Moreover Nisini et al. 
(1995) suggested that TTS have higher Mwind (based on 
TTS data from Giovanardi et al. 1991) than one would 
expect from extrapolating the M^ind vs. Lboi relationship 
of the HAEBES to lower luminosities. Here Lboi is the 
bolometric luminosity of the source and thus includes not 
only the stellar but the disk and any halo component as 
well. Putting aside, for the moment, the claim by Calvet 
(1997) and others that the hydrogen emission lines in any 
event do not arise in the wind, the suggestion that TTS are 
somehow more "efficient" at losing mass than HAEBES is 
difficult to reconcile with our findings: we see a smooth 
transition in the relationship between wind strength and 
accretion across the TTS/HAEBES boundary (see Fig. 5). 
Moreover we also find a smooth relationship, between Ha 
line strength and infrared excess. It appears, however, that 
the TTS mass-loss rates calculated by Giovanardi et al. 
(1991) are in any event highly unreliable (Natta, private 
communication) . In fact many of the TTS listed by Giova- 
nardi et al. (1991) also had their mass loss rates estimated 
by Hartigan et al. (1995). These authors found values typ- 
ically 10 times lower than Giovanardi et al. (1991), imply- 
ing the "efficiency" of mass loss is not higher amongst TTS 
than HAEBES and implying a much smoother transition 
in Mwind vs. Lboi than found by Nisini et al. (1995; ibid 
Fig. 4). It is also worth remarking that a direct compar- 
ison of the Lire and L» values in Nisini et al. (1995) and 
ours reveals some distinct differences. It appears, however, 
that the calculated values of L, in Nisini et al. (1995) did 
not take account of the appropriate bolometric correction 
to Mv 

As can be seen from Figs. 4 & 5, there is a tendency for 
L[oi], Lho) and Lire to scale with the stellar luminosity, L«. 
In particular Fig. 4a shows that there is a good correlation 
between Lire and L* and it is this relationship which, to 
some degree, masks whether L[oi] and Lho depend fun- 
damentally on Lire or L*. The form of the dependence is 
critical to understanding the origin of the forbidden line 
emission in HAEBES: Catala (1989) proposed a chromo- 
spheric origin for the winds from HAEBES that superfi- 
cially might be in keeping with the dependence of L[oi] , for 
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example, on L*. According to that paradigm, the observed 

correlations with Lire would then be coincidental and due 
to the tendency for the total amount of circumstellar dust 
to scale with L*. Such an explanation, however, does not 
explain the tighter correlation we observe between L[oi] 
and Lire than between Ljqi] and L». Moreover, a chromo- 
spheric model seems incapable of explaining why a num- 
ber of HAEBES show infrared excess luminosities that arc 
much greater than the underlying stars' bolometric lumi- 
nosities (Hamann & Persson 1992). Accretion must play 
a fundamental role as it is also impossible to imagine how 
such a situation can arise in a purely passive reprocessing 
disk or envelope. In fact, it may even be argued that the 
"reprocessing fraction" , particularly for the earliest stars, 
is lower amongst HAEBES than TTS due to dust clearing 
(Hillenbrand et al. 1992; Hamann & Persson 1992). 

A final point should be made in connection with the 
apparent dependence of wind properties on L* amongst 
HAEBES. In the case of the TTS, Cabrit et al. (1990) 
found a poor correlation of Lpj] and Ljja with TTS pho- 
tospheric luminosity in contrast to what is found for the 
HAEBES. The probable explanation for this effect is that 
in the TTS sample one is dealing with a large number of 
stars within a rather narrow luminosity (and mass) range. 
The sample of Cabrit et al. (1990) ranges in spectral type 
from Ml to G2 or a mass range of roughly 0.4 - 1 Mq. The 
scatter in such a sample of Lire, for a given narrow range 
in L*, is then large enough that the stronger dependence 
on Lire of L[oi] , for example, is readily seen. 

It is interesting to compare the frequency of the occur- 
rence of molecular (CO) outflows with the near-infrared 
colour, or effectively as one can see from Fig. 1, with Hil- 
lenbrand group. While none of our Group III stars are ob- 
served to be accompanied by molecular outflows and only 
2 out of our 21 Group I stars arc, 6 out of our 9 Group II 
stars are found to be associated with molecular outflows. 
Note that we used the molecular outflow database of Fukui 
et al. (1993) to check for associations. There thus appears 
to be a direct link between degree of embeddedness of the 
star and the presence or absence of large scale outflow 
phenomena. The Group II stars' association with molecu- 
lar outflows, and in some cases optical jets (Mundt & Ray 
1994), may be due to their relatively early evolutionary 
state, as they disperse the dusty cocoons from which they 
formed. Strom et al. (1993) have suggested that the ces- 
sation of the outflow phase in TTS is intimately linked to 
the circumstellar disk surrounding the star becoming opti- 
cally thin. The Group III HAEBES may possess optically 
thin circumstellar material and colder dust at large radii, 
a remnant of an earlier phase, responsible for the longer 
wavelength IR excesses observed in such stars. 

4-3. Correlations with the circumstellar mass 

In Fig. 6a we showed a plot of the luminosity of the 
[OIJA6300 line against the calculated circumstellar mass. 



There is a clear correlation between the two quantities as 

is also the case for the infrared excess (Lire) and the cir- 
cumstellar mass (see Fig. 6b and also Table 4). Obviously 
to an extent such correlations reflect the fact that the cir- 
cumstellar mass tends to scale with the stellar luminosity 
L* and, given that Lpij and Lire depend on L«, are to 
be expected. Nevertheless, Table 4 shows that the corre- 
lations of L[oi] and Lire with Mes arc stronger than with 
L« suggesting some dependence of these quantities on Mes 
alone. 

The origin of the mm emission in HAEBES is not clear 
and this is partly due to the poor resolution of single dish 
mm measurements such as those reported here. As was 
stated in the Introduction, the SEDs of HAEBES have var- 
iously been modeled by extended dusty envelopes without 
disks (Berrifli et al. 1992; Di Francesco et al. 1994; Mirosh- 
nichcnko et al. 1997) as well as by disk plus envelope com- 
binations (Hillenbrand et al. 1992; Natta 1993). In order to 
determine whether there is a disk contribution to the mil- 
limeter emission, higher resolution mm-array observations 
are needed. To date, however, these have produced con- 
flicting results (e.g., Di Francesco et al. 1997; Mannings & 
Sargent 1997) In any event it is clear that if a major frac- 
tion of the circumstellar material inferred from our mm 
observations is conflned to scales of typically a few arcsec- 
onds then, given its mass (see Table 5), one would expect 
much higher visual extinctions towards these sources than 
observed where the material to be distributed spherically. 
We emphasize, however, that even if disks are present, as 
we believe they are, this does not exclude the possibility 
that some of the mm emission almost certainly arises from 
an extended dusty halo (Di Francesco et al. 1994). 

The correlations we observe in Fig. 6 do not tell us any- 
thing about the origin of the mm emission. They would 
be expected in the event that most of the mm emission 
came from an accretion disk but equally they might be an- 
ticipated if the mm emission largely derived from a dusty 
envelope that in turn fed a disk. 

Finally we add a note of caution that, given our beam 
size, our mm flux measurements may include contribu- 
tions from unrelated dust, and/or embedded companions, 
at large distances from the HAEBES. It is interesting 
to note, however, that our circumstellar mass estimates 
are close to those derived by Mannings & Sargent (1997), 
from their much higher resolution observations, in at least 
the two stars common to both samples (i.e. AB Aur and 
HD 163296). 

5. Conclusions 

The work of Cohen et al. (1989), Cabrit et al. (1990), Ed- 
wards et al. (1993), Hartigan et al. (1995) and others has 
demonstrated a strong link between accretion and mass- 
loss in the case of the CTTS. Using various indicators of 
Mace, the mass accretion rate, and Mwind, the mass- loss 
rate, they have shown that the mass loss rate scales with 
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the accretion rate and that the mass-loss most likely de- 
rives from an accretion-driven disk flow. The motivation 
for the present, and related (Corcoran & Ray 1997a), work 
is to test whether the same basic inflow/outflow process 
occurs amongst the HAEBES. The approach in this pa- 
per has been to test whether the same relationships seen 
in the TTS between indicators of mass loss and mass ac- 
cretion carry over to their higher mass counterparts. Ob- 
viously if clear parallels can be drawn between the two 
groups this would not only bolster the idea that a com- 
mon inflow/outflow mechanism operates amongst the TTS 
and HAEBES but indirectly, that disks, as in the case of 
CTTS, surround HAEBES (or at least those with sub- 
stantial infrared/mm excesses). Certainly we expect that 
the disk paradigm can be extended to some intermediate 
mass young stars given the presence of dusty disks/rings 
around main sequence A-type stars like /? Pic (Burrows et 
al. 1997). The crucial question is rather how common are 
disks around HAEBES? The data and analysis we have 
presented here shows that in the case of the HAEBES: 

— Relative strengths of [OI]A6300 and Ha lines, as deter- 
mined from absolute equivalent widths, typically scale 
with reddened corrected ncar-IR colours, H-K and K- 
L for the Hillenbrand et al. (1992) Group I and II 
HAEBES. As discussed by Edwards et al. (1993) for 
the CTTS, these near-IR coloins sample the tempera- 
ture distribution of the inner system dust around the 
star. As the HAEBES near-IR colours increase above 
the values expected from a passive reprocessing or ge- 
ometrically thin accretion disk, it is found that the 
relative wind strength increases in an exactly analo- 
gous manner to the CTTS. Since Hartigan et al. (1995) 
have shown that, at least in the case of the CTTS, the 
reddening corrected IR colours correlate well with the 
degree of veiling and hence indirectly with the relative 
accretion rate, this strongly suggests that the same 
basic inflow/outflow process occurs in both CTTS and 
HAEBES. The idea that relative accretion is related 
to near-IR colour in HAEBES is further reinforced by 
the finding that for these stars K-L is correlated with 
Lire/L, where Lire is the infrared excess. Lire, if one 
ignores the effects of reprocessing, should be a rough 
measure of accretion. 

— There are indications of a small population of Group 
I stars with no detected [OIJA6300 emission and the 
levels of which are well below that expected on the 
basis of their near-IR colours. It seems unlikely that 
their infrared excesses can be explained by the presence 
of companions and the absence of [OI]A6300 emission 
may simply be due to variability as is certainly true 
in a number of cases. Another possibility is that their 
disks are currently quiescent and merely reprocessing 
starlight. 

— All Group II HAEBES have near-IR colours that 
are much redder than those of a passive reprocess- 
ing or geometrically thin accretion disk. As noted 



by Hillenbrand et al. (1992) their overall spectral 

energy distribution is best explained in terms of a 
star-l-disk-l-envelope model. These are found to be the 
the stars with the strongest evidence for accretion. 
We note that the Group II stars are, proportionately, 
by far the most likely sources of extended optical jets 
amongst all HAEBES (Corcoran & Ray 1997a). 

— No forbidden line emission was detected from any of 
our Hillenbrand et al. (1992) Group III stars in keep- 
ing with the small values of their near-infrared colours 
(i.e. close to that expected from a main sequence star 
of the same spectral type). Moreover the equivalent 
width of their Ha emission is < IsA in all cases. 
We propose that their lack of forbidden line activ- 
ity and, in particular such an Ha equivalent width 
limit, may serve as a useful criteria when searching for 
members of the Group HI HAEBES sub-class. That 
said, we should warn the reader that the weak-line 
HAEBES, "WHAEBES", defined in this way are not 
exactly identical to the Group III stars as defined 
by Hillenbrand et al. (1992). The latter are charac- 
terised in terms of the near-infrared slope of their 
SEDs whereas the WHAEBES will also include those 
Group I /II stars with low line activity (e.g. BF Ori, 
MWC 137, HD 150193), the status of which is still 
unclear. While it seems some Group III stars are not 
pre-main sequence stars, nevertheless it appears likely 
that a number represent the HAEBES analogues of 
the weak-line T Tauri stars. Interestingly Group HI 
stars are predominantly early B-type stars. Such stars 
evolve rapidly and will therefore not have wandered 
far from their formation sitc;s. On the other hand the 
Group I and II stars are predominately of later spectral 
type. This begs the question where are the Group HI 
late-type HAEBES? Such stars should be searched for, 
perhaps many of them lie well beyond the boundaries 
of their parent molecular cloud. 

— Using the luminosity of the [OIJA6300 emission line as 
an indirect measure of Mwind and the IR excess lumi- 
nosity over 0.7-10.2/im as a measure of Mace, a strong 
correlation is found between accretion and mass-loss. 
While the [OIJA6300 line luminosity also tends to in- 
crease with stellar luminosity, as one might expect, the 
correlation with IR excess luminosity is stronger im- 
plying that accretion rate is an important factor in 
determining the strength of the [01] emission. 

— Plotted to the same axes as our HAEBES data, the 
L[oj] vs. Lire of the sample of 36 TTS from Cabrit 
et al. (1990) extends the correlation between the two 
variables over 5 orders of magnitude in luminosity. The 
smooth progression from PMS stars of less than a solar 
mass all the way to B-type stars with masses ~ IOMq 
or more provides convincing evidence that the origin 
of the Mwind « Mace relationship is the same process 
in both groups of stars. Given the compelling evidence 
in the case of CTTS that matter is accreted via a disk. 
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one must infer that accretion disks also surround many 
HAEBES. 

— Millimeter continuum observations of HAEBES reveal 
circumstellar masses typically in the range 0.02-1.0 
M0, assuming a gas to dust ratio of 100:1 and opti- 
cally thin emission. It is likely in a number of cases 
that these masses are somewhat overestimated due to 
our large beam size (30"). If a significant proportion 
of the mm flux comes from a compact region on scales 
less than a few arcseconds then the masses inferred 
would imply much higher optical extinction towards 
these stars than observed unless the dust is not dis- 
tributed spherically. Recent observations support the 
idea that the emission is compact (Mannings & Sar- 
gent 1997) and this again suggests that disks may be 
present around many HAEBES. The amount of cir- 
cumstellar mass, Mcs, present is slightly better corre- 
lated with the luminosity of the forbidden line emis- 
sion than with the photospheric luminosity of the star. 
Again this indicates that a disk or a disk-|-halo com- 
bination may be fueling HAEBES with forbidden line 
emission. 
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